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Quasi-crystals formed by charged mesoscopic dust grains (dust lattices), observed since hardly a 
decade ago, are an exciting paradigm of a nonlinear chain. In laboratory discharge experiments, 
these quasi-lattices are formed spontaneously in the sheath region near a negative electrode, usually 
at a levitated horizontal equilibrium configuration where gravity is balanced by an electric field. It 
is long known (and experimentally confirmed) that dust-lattices support linear oscillations, in the 
longitudinal (acoustic mode) as well as in the transverse, in plane (acoustic-) or off-plane (optic-like 
mode) directions. Either due to the (typically Yukawa type) electrostatic inter-grain interaction 
forces or to the (intrinsically nonlinear) sheath environment, nonlinearity is expected to play an 
important role in the dynamics of these lattices. Furthermore, the coupling between the different 
modes may induce coupled nonlinear modes. Despite this evidence, the elucidation of the nonlinear 
mechanisms governing dust crystals is in a rather preliminary stage. In this study, we derive a set 
of (coupled) discrete equations of motion for longitudinal and transverse (out-of-plane) motion in a 
one dimensional model chain of charged dust grains. In a continuum approximation, i.e. assuming 
a variation scale which is larger than the lattice constant, one obtains a set of coupled modified 
Boussinesq-like equations. Different nonlinear solutions of the coupled system are discussed, based 
on localized travelling wave ansatze and on coupled equations for the envelopes of co-propagating 
quasi-linear waves. 
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I. INTRODUCTION 
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Recent studies of various collective processes in dust contaminated plasmas (DP) have been of significant interest 
in relation with linear and nonlinear waves which are observed in laboratory and space plasmas. An issue of particular 
7-H , importance is the formation of strongly coupled DP crystals by highly charged dust grains, for instance in the sheath 
region above a horizontal negatively biased electrode in experiments 0, . Low- frequency oscillations may occur in 
these mesoscopic dust grain quasi-lattices, in both longitudinal (acoustic mode) Q and transverse (in-plane shear 
acoustic mode, off-plane optic-like mode) directions, as theoretically predicted and experimentally observed (see in 
Ref. l] for a review). 

In this paper, we focus on the nonlinear description of dust grain displacements in a one-dimensional dust crystal, 
which is suspended in a levitated horizontal equilibrium position where gravity and electric (or, possibly magnetic 0) 
forces balance each other. Considering the coupling between the horizontal (~ x) and vertical (off-plane, ~ z) degrees 
of freedom, and an arbitrary inter-grain interaction potential form U(r) (e.g. Debye or else) and sheath potential 
O \ *&{ z ) ( n °t necessary parabolic), we aim in deriving a set of equations which should serve as a basis for forthcoming 
studies of the nonlinear behaviour of longitudinal and transverse dust lattice waves (LDLWs, TDLWs) propagating 
in these crystals. The relation to recent studies of a similar scope (here recovered as special cases) is also discussed. 

Or 
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II. THE MODEL 

Let us consider a layer of charged dust grains (mass M and charge q, both assumed constant for simplicity) of 
lattice constant r . The Hamiltonian of such a chain is of the form 

where r„ is the position vector of the n— th grain; U nm (r nm ) = q<fi(x) is a binary interaction potential function related 
to the electrostatic potential <fi(x) around the m— th grain, and r nm = r„ — r m \ is the distance between the n— th and 
m — th grains. The external potential $ ext (r) accounts for the external force fields in which the crystal is embedded; 
in specific, <& ex t takes into account the forces acting on the grains (and balancing each other at equilibrium, ensuring 
stability) in the vertical direction (i.e. gravity, electric and/or magnetic forces); it may also include the parabolic 
horizontal confinement potential imposed in experiments for stability Q as well as, for completeness, the initial laser 
excitation triggering the oscillations in experiments. 




FIG. 1: Dust grain vibrations in the longitudinal (~ x) and transverse (~ z) directions, in a Id dust lattice. 



A. 2d equation of motion 

Considering the motion of the n— th dust grain in both the longitudinal (horizontal, ~ x) and the transverse (vertical, 
off-plane, ~ z) directions (i.e. suppressing the transverse in-plane - shear - component, ^ SO that T n — {x n ^ 
we have the two-dimensional (in x, z) equation of motion 

\~a¥~ ~d±) = dr F ^t( r «) = «E(r„) + F ext (r n ) , (l) 

where Ej(x) = —d<j)(r)/dxj is the (interaction) electrostatic field and F ex tj — —d$ e xt(x)/dxj accounts for all external 
forces in the j— direction (j = 1/2 for Xj = x/z); the usual ad hoc damping term was introduced in the left-hand-side 
of Eq. involving the damping rate v due to dust -neutral collisions. 

B. Nonlinear vertical confining potential 

We shall assume a smooth, continuous variation of the (generally inhomogeneous) field intensities E and/or B, as 
well as the grain charge q (which may vary due to charging processes) near the equilibrium position z — 0. Thus, we 
may develop 

E(z)^E + E' z+ X -E'^z 2 + 
B(z)^B + B' z+ ^B{z 2 + 
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and 



q(z) « q + q' z + -q'^z 2 + 

where the prime denotes differentiation with respect to z and the subscript '0' denotes evaluation at z — Zq, viz. 
Eq = E(z = zq), E' = dE(z) /dz\ z=Zo and so forth. Accordingly, the electric force F e = q(z)E(z) and the magnetic 
force F m — —d(mB)/dz = —2aBdB/dz (where the grain magnetic moment fi is related to the grain radius a and 
permeability /i via m = (fi — l)a 3 B/(ji + 2) = aB |(|), which are now expressed as 

1 



F e (z) « q E + {q E' + q' a E Q ) z + ^{q E'J + 2q' E' a + q{E ) z 2 + 



and 



F m (z) « -2aB Q B' - 2a(B' Q 2 + B Q B" Q ) z - a(B B%' + 3B' B" ) z 2 
may be combined to give 

F + F — 

oz 

where we have introduced the phenomenological potential $(z) 
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The definitions of = (9 J $(z)/9z : '| z _ z = — (qEo)^ ' + a^B 2 )^ (here, the superscript within parenthesis 
obviously denotes the order in partial differentiation; j = 1, 2, ...) are obvious: 

= ~{qE) + a{B 2 y = -q E + 2aB B' 



$(2) = 



4> 



(3) 



-(qE y + a(B 2 )% 

-(q' E + q E' Q ) + 2a{B'l + B Q B'J) 

-(q%E + 2q' E' + q E{) + 2a(3B' B% + B Q B'J') , 



(3) 



and so forth. Obviously, & e xt = $ — Mgz. The (vertical) force balance equation d$ ex t/dz = 0, viz. 

Mg = q E - 2aB B' , 

is satisfied at equilibrium. 



III. DISCRETE EQUATIONS OF MOTION 

Assuming small displacements from equilibrium, one may Taylor expand the interaction potential energy U(r) 
around the equilibrium inter-grain distance Itq = \n — m\ro (between l—th order neighbors, I = 1,2, ...), i.e. around 
5x n Ri and 5z n w 0, viz. 

(x n X m ) 3 

r— I |n— m\ro 

where I' denotes the degree of nonlinearity involved in its contribution: I' = 1 is the linear interaction term, V = 2 
stands for the quadratic potential nonlinearity, and so forth. Notice that the inter-grain distance r = [(x n — x m ) 2 + 
(z n — Zm) 2 ] 1 / 2 also needs to be expanded, i.e. near \x n — x m \ = Iro and z n — z m = 0, so that dll(r)/dxj = 

(dU(r)/dr)(dr/dxj) ~ .... Obviously, 5x n = x n — x^ and 6z n — z n — Zn denotes the displacement of the n— th 



Id' 

u(r nm ) = J2 jrr 

ii n 
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grain from the equilibrium position (xn \ Zn^) = (wo, 0). Retaining only nearest-neighbor interactions (I = 1), we 
obtain the coupled equations of motion 



d 2 (8x n ) d(Sx n ) 



dt 2 



dt 



u>l L (fen+l + Sx n -i - 2Sx n ) 



O30 



(Sx n+ i - 5x n ) 3 - (8x n - SXn^xY 



ao2 



(SXn+1 - Sx n ) 2 - (Sx n - 5x n -i) 2 

(Sz n+1 - Sz n ) 2 - (Sz n - 5z n -i) 2 



-a 12 



(5x n+1 - 5x n )(Sz n+1 - bz n ) 2 - (Sx„ - 5x n -x){8z n - 8z n -if 



and 



d 2 (Sz n ) d(Sz n ) 2 2 

+ v — - — = w 0T (2dz n - oz n+1 + dZn-i) - lo dz n 



dt 2 



dt 



- K x (5z n y - K 2 {8z n y + 



3 . a 02 



ro 



(8z n+1 - Sz n ) 3 - (Sz n - 6z n -iY 



2a 2 



(Sx n+ i - Sx n )(Sz n+1 - Sz n ) - (Sx n - 6x n -i)(6z n - 5z n -{) 



a\2 



(Sx n+ i - Sx n ) 2 (Sz n+1 - Sz n ) - (Sx n - fa„„i) 2 ((5z„ - 6z n -i) 



where we have defined the longitudinal/transverse oscillation characteristic frequencies 

"l L = U"(r )/M , Lu 2 tT = - U'(r )/(Mr ) , 
(both assumed to be positive for any given form of interaction potential U) and the quantities 

«20 = ~ U"'(r ) , a Q2 = --^ [U'(r ) - r U"{r )] , 
«30 = 7^7 U""(r ) , a 12 = - J3 [U'(r ) ~ r U"(r ) + r 2 \ U"'(r )] , 



6M 



Mrl 



(4) 



(5) 



(6) 



(7) 



which are related to coupling nonlinearities. The gap frequency to g and the nonlincarity coefficients K\ and K 2 are 
related to the form of the sheath environment (i.e. the potential $) via 



lu 2 = $ (2) /M, X 1 = $ (3) /(2M), K 2 = $ (4) /(6M) . 



(8) 



Obviously, the prime denotes differentiation, viz. U"(ro) — d 2 U(r) / dr 2 \ r _ r and so on. In the above equations of 
motion, we have distinguished the linear contributions of the first neighbors from the nonlinear ones, i.e. the first 
line in the right-hand-side from the remaining ones, in both equations. Note that all of the coefficients are defined in 
such a way that they bear positive values for Debye-type interactions, i.e. if Ud(t) — (q 2 /r) exp(— r/\]j) (Xd is the 
effective Dcbye length) since odd/even derivatives are then negative/positive; however, the sign of these coefficients is 
not a priori prescribed for a different interaction potential U(r). Indeed, we insist on expressing all formulae in such a 
manner that a different interaction law may easily be assumed in a "plug-in" manner; in particular, even though the 
Debye potential Ud is widely accepted in DP crystal models, we think of the modification of U when one takes into 
account a magnetic field £| or the ion flow towards the negative electrode 0- Nevertheless, we provide the explicit 
form of the coefficients a.y defined above for a Debye potential, for clarity, in the Appendix. 

Upon careful inspection of the discrete equations of motion above, one notices that the lowest order nonlinearity in 
the longitudinal motion is due to the intcrgrain interaction law, while nonlinearity in the vertical motion is primarily 
induced by the coupling to the horizontal component (and, to less extent, by interactions). 
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IV. CONTINUUM APPROXIMATION 



Adopting the standard continuum approximation, we may assume that only small displacement variations occur 
between neighboring sites, i.e. 

du 1 2 9 2 u 1 3 d 3 u 1 4 d 4 u 
Sx n±1 ^u±r - + -r Q —± -r — + -r — ± 

where the (horizontal) displacement 5x n (t) is now expressed via a continuous function u = u(x,t). The analogous 
continuous function w = w(x,t) is dchncd for 8z n (t). 

One may now proceed by inserting this ansatz in the discrete equations of motion and carefully evaluating 

the contribution of each term. The calculation, quite tedious yet perfectly straightforward, leads to a set of coupled 
continuum equations of motion in the form 

c? 



u + v u - c| u xx - ^1 rjj u xxxx = - 2 a 2 a Tq u x u xx + 2 a 02 w x w xx 
A \t... ^2 „, i o„„ ... .. 1 i q „4 / ^2 , 



- 012^0 [(wj) u xx + 2w x w xx u x ] + 3a 30 r (u x ) u xx , (9) 



U> + 1/ W + c| W xx + j| r w xxxx + UJg w = - K 1 w 2 - K 2 w 3 

+ 2 a 02 rl (u x w xx + w x u xx ) 
+ 3 ao2 Tq (w x ) w xx — ai2 r o [(u x ) w xx + 2u x u xx w x ] , (10) 

where higher-order nonlinear terms were omitted. We have defined the characteristic velocities cl — wo.l and 
ct = Wo,T r o! the subscript x denotes partial differentiation, so that u x u xx = (u x ) x /2 and (u x ) 2 u xx = (u x ) x /3. 
Remember that the gap frequency to g and the coefficients K\ and Ki are related to the form of the sheath electric 
and/or magnetic potential via JSJ above, viz. F e i — M g — M uig z — K\ z 2 — z 3 . 



V. RELATION TO PREVIOUS RESULTS - DISCUSSION 



As a matter of fact, all known older results are based on equations which are readily recovered, as special cases, 
from Eqs. JH) and (JSJ and/or their continuum counterparts © and (JTUJ. In particular, the coupled Eqs. (1) and 
(2) in Ref. y| are exactly recovered from and J5J), upon neglecting 030, ai2, K\ and K2 and then evaluating all 
coefficients for a Debye-type potential. 

Upon switching off the coupling (i.e. setting w — > 0), Eq. © above recovers exactly the nonlinear Eq. (13) in Ref. 
|lf)j |. which was therein shown to model (nonlinear) longitudinal dust grain motion in terms of (either Korteweg-de 
Vries- 0, or Boussinesq-type) solitons; also see Eq. (2) in ^l| (treating the formation of asymmetric envelope 
modulated LDLWs) and Eq. (2) in ^3] (keep only first-neighbor interactions therein, to compare). In a similar 
manner, considering purely transverse motion (i.e. setting u — > 0) Eqs. (J5J and (|l(Jf> herein recover exactly the 
nonlinear Eqs. (7) and (8) in Ref. 01, where they were shown to model the amplitude modulation of TDLWs which 
is due to the sheath nonlinearity. Finally, needless to say, the linear limit recovers exactly the known equations of 
motion for cither purely longitudinal or purely transverse motion (i.e. considering ay = Kj =0, Vi,j). 

An exact treatment of the coupled evolution Eqs. - or, at least, the continuum system 11U[) - seems quit e 

a complex task to accomplish. Even though Eq. I© may straightforward be seen as a Boussinesq-type equation [l(| , 
which is now modified by the coupling, its transverse counterpart (|10|l (for u — > 0, say) substantially differs from 
any known nonlinear model equation, bearing known exact solutions. Therefore, we shall limit ourselves to reporting 
this system of evolution equations, for the first time, thus keeping a more thorough investigation (analytical and/or 
numerical) of their fully nonlinear regime for a later report. 



VI. COUPLED-MODE MODULATED WAVE PACKETS 



In order to gain some insight regarding the influence of the mode-coupling on the nonlinear profile of the dust 
lattice waves, we may consider the effects which come into play when the amplitude of the LDLWs and the TDLWs - 
which are initially uncoupled in the small amplitude (linear) limit - is increased to a slightly finite (i.e. non negligible) 
value, thus allowing for a weak coupling between the two modes and a tractable appearance of the signature of the 
(weak) nonlinearity in the dynamics. 
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The standard way for such an approach is via the introduction of multiple space and time scales, viz. Xq, Xi,X 2 , ... 
and T , T l7 T 2 , where X n = e n x and T n — e n t (e <C 1 is a smallness parameter). The solutions are expanded as: 
u = eui + e 2 u 2 + ... (plus an analogous expression for w). The technical details of the calculation are described e.g. 
in p"H and will be omitted here. We shall apply this reductive perturbation technique to the system obtained from 
Eqs. 1 1 0|t by keeping only the lowest-order nonlinear terms (i.e. omitting the last line in both equations); we set 
po = 2a2o^o an d = 2ao2?"o f° r simplicity. Note the inevitable (and qualitatively expected) complication of the 



le (a 
14]- 



calculation due to the different dispersion laws in the two modes 

The first-order (~ e) equations are uncoupled and may be solved by assuming {ui,wi} — {ip^jip^} 

+[{ipL, "0t} e~xpi(kx — cot) + c.c] (complex conjugate). Upon substitution, we obtain ip^ = 0; the remaining (3) 
amplitudes are left arbitrary. This readily yields the known dispersion relations 

uji+iuu L = c? L e(i-^^j , j r + lV u T = ui~? T k 2 (\- ] ^^, (ii) 

for the (acoustic) LDL and the (optical-like) TDL mode respectively. 

The 2nd-order (~ e 2 ) equations contain secular (lst-harmonic forcing) terms, whose elimination imposes a pair of 
conditions in the form: j/dT\ + v g j d^j/dXi = (where j £ {1, 2} = {L,T} in the following), implying that the 
amplitudes ^ j travel at the (different) group velocities v g j = dujj(k)/dk. See that v g ^r — wr'(^) < (the TDLW is 
a backward wave), as immediately obtained from i|llb). The remaining system is then solved for the 0th and the 2nd 
harmonic amplitudes (in e 2 ) 14]; the solution finally obtained is of the form: 

Sx n (t) « u(x, t) « e[-0o + -01 expi(kx — lo\1) + c.c] + e 2 u 2 exp 2i(kx — LO\t) + c.c] + C(e 3 ) 

Sz n (t) « w{x,t) w e[if} 2 expi(kx - u) 2 t) + c.c] + e 2 {w { 2 0) + [wf* exp2i(kx ~ w 2 t) + c.c.]} + C(e 3 ) . (12) 

We henceforth denote the significant amplitudes and by ^o, and ^ 2 respectively. The 2nd order 

correction amplitudes are 

u ( 2) = tk3 Po n -hon m = w = _ j_ ( KM + 2ihok ^ 2 ) , (is) 

where = ~ c | r o^ 4 + 2iwi and D% = ~ 3cj 2 + c^rpfc 4 + 2wlot- The contributions u 2 x \ w 2 an d "2 are left 
arbitrary by the algebra and were thus set to zero. 

Proceeding to the 3rd-order (~ e 3 ) equations, the elimination of the secular terms together with zeroth order 
equations provide three explicit conditions, for ^0,1,2- After some tedious algebra, these take the form 

i {^ +Vg ' 1 ^£) + Pl ^XJ + l*i| 2 *i + Qi 2 |* 2 | 2 *i +(Qo,i*i + Qo,2* 2 ) H^ + ^i =0 

(9*9 5*9 \ 9 2 *9 

-q^T + J + P2 -g^i- + Q22 |* 2 | 2 *2 + Q21 |*l| 2 *2 + H-2 = 

«r - c|)|^ = - Po fc 2 |*i| 2 + h a k 2 \V 2 \ 2 + C , (14) 

where C is an integration constant (to be determined by the boundary conditions). The linear dispersion terms Pj 
are related to the (curvature of) the dispersion relations ifTTjl as Pj = ujj"{k) {j = 1,2); the group velocities v g j 
were defined above [l^. The nonlinearity coefficients Qij (i = 0,1,2, j = 1,2) and the 'peculiar' contributions Hj 
(involving cross-terms in 'J 2 **) are too lengthy to report here Observe that, once C is determined, one may cast 
Eqs. I|14|l into the form of a (modified, asymmetric) system of coupled nonlinear Schrodinger equations (CNLSE). 
Note that we have avoided the usual envelope (Galilean) transformation {x,t} — > {x ~ v g ^t,t}, since it does not 
simplify this (asymmetric, with respect to 1 <-» 2) system. Finally, let us point out, for rigor, that the results in 
[ill ] and 0] are exactly recovered, from both l|13|) and l|14|) . in the appropriate - uncoupled mode - limits (namely, 
* 2 — > and *i — > 0, respectively, for LDLWs and TDLWs). 

Despite the obvious analytical complication, the physical mechanism underlying the above results is rather trans- 
parent. There is an energy pumping effect between the zeroth-harmonic longitudinal (displacement) mode *0j first 
put forward in (for LDLWs, yet long known in solid state physics 01) anc l th e modulated (low- frequency) LDL 
and (high-frequency) TDL mode(s). Note the strong misfit (asymmetry) between the dispersion laws dominating 
the coupled modes, despite which - regretfully - no simplifying assumption may be analytically carried out in this 
continuum model. 
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VII. CONCLUSION 



We have put forward a comprehensive nonlinear model for coupled longitudinal-to-transverse displacements in a 
horizontal dust mono-layer, levitated in a sheath under the influence of gravity and an electric and/or magnetic field. 
All of the above results are generic, i.e. valid for any assumed form of the inter-grain interaction potential U(r) and 
the sheath potential $, and will hopefully contribute to the elucidation of the grain oscillatory dynamics in dust 
crystals. 

Appendix: Form of the coefficients for the Debye interaction potential 

Consider the Debye potential (energy) Uo(r) — q(t>D( r ) = Q 2 e~ r l Xo jr. Defining the (positive real) lattice parameter 
k = tq/Xo, one straightforward has 



U D {r ) = -Ti-e — — > U D (r ) = +- r e 



2 



U D { r Q> ~ Af, K 4 U D VO) ~ + \$ d e K 5 

where the prime denotes differentiation and I = 1, 2, 3, ... is a positive integer. Now, combining with definitions JfjJE}, 
we have: 

.2 2<7 2 K 1 + K + K 2 /2 _ 2 2 2 2 gr 2 1 + K - „ „ 



p = 2a 20 K 3 \ 3 D = e- K f i + 1 + J + ) , h = 2a m n 3 X 3 D = e~ K ( - + 1 + J 



MX D \k 2 6 J u MX D \k 3 

a-30 = q r* e_K A f« 4 + 4k3 + 12 ^ 2 + 24k + 24 ^1 > fl i2 = — tttkt er K \ ( n 3 + 5k 2 + 12k + 12 



6MX 5 D k 5 V J ' 2M A|, k 5 



Of course, all known previous definitions of (some of) these coefficients (for nearest neighbour interactions; see in the 
references cited in the text) are exactly recovered. Note, finally, that n is of the order of (or slightly higher than) unity 
in experiments; therefore, all coefficients turn out to be of similar order of magnitude, as one may check numerically. 



Acknowledgments 

This work was supported by the SFB591 (Sonderforschungsbereich) - Universelles Verhalten gleichgewichtsferner 
Plasmen: Heizung, Transport und Strukturbildung German government Programme. Support by the European Com- 
mission (Brussels) through the Human Potential Research and Training Network via the project entitled: "Complex 
Plasmas: The Science of Laboratory Colloidal Plasmas and Mesospheric Charged Aerosols" (Contract No. HPRN- 
CT-2000-00140) is also acknowledged. 



[f] P. K. Shukla and A. A. Mamun, Introduction to Dusty Plasma Physics (Institute of Physics, Bristol, 2002). 

[2] G. E. Morfill, H. M. Thomas and M. Zuzic, in Advances in Dusty Plasma Physics, Eds. P. K. Shukla, D. A. Mendis and 

T. Desai (Singapore, World Scientific) p. 99. 
[3] F. Melands0, Phys. Plasmas 3, 3890 (1996). 

[4] V. V. Yaroshenko, G. E. Morfill and D. Samsonov, Phys. Rev. E 69, 0I64I0 (2004). 

[5] D. Samsonov, A. V. Ivlev, R. A. Quinn, G. Morfill and S. Zhdanov, Phys. Rev. Lett. 88, 095004 (2002). 
[6] J. D. Jackson, Classical Electrodynamics (John Wiley and Sons, New York, 1963). 

[7] A. M. Ignatov, Plasma Physics Reports 29, 296 (2003); I. Kourakis and P. K. Shukla, Phys. Lett. A 317, 156 (2003). 



8 



[8] A. V. Ivlev, S. K. Zhdanov, and G. E. Morfill Phys. Rev. E 68, 066402 (2003). 
[9] L. Stenflo, N.L. Tsitsadze and T. D. Buadze, Phys. Lett. A 135 (1), 37 (1989). 
[10] I. Kourakis and P. K. Shukla, Eur. Phys. J. D, in press (2004) (DOI: 10.1140/epjd/e2004-00025-3). 

[11] I. Kourakis and P. K. Shukla, Modulated wave-packets associated with longitudinal dust grain oscillations in a dusty plasma 

crystal, Phys. Plasmas, 11 (5) (in press, sch. May 2004). 
[12] K. Avinash, P. Zhu, V. Nosenko and J. Goree, Phys. Rev. E 68, 046402 (2003). 

[13] I. Kourakis and P. K. Shukla, Weakly nonlinear vertical dust grain oscillations in dusty plasma crystals in the presence of 

a magnetic field, Phys. Plasmas, submitted (2004). 
[14] The tedious details are left to be reported in a lengthier report, in preparation. 

[15] See that, once damping is taken into account (via v 7^ 0), all dispersion-related quantities - namely u)j, v Bt j and Pj 
(j = 1,2 = L,T) here - bear an imaginary part. All of the corresponding expressions - readily obtained from Eq. 1111 
and thus omitted here - were exactly recovered by the (tedious) algebra within our perturbative scheme; recall, in passing, 
that ^o t i t 2 are complex numbers, by definition. 

[16] I. Kourakis and P. K. Shukla, http://arxiv.org/abs/cond-mat/0402007, to appear in Phys. Plasmas 11, (2004). 

[17] A. Tsurui, Progr. Theor. Phys. 48, 1196 (1972); N. Flytzanis, St. Pnevmatikos and M. Remoissenet, J. Phys. C: Solid 
State Phys. 18, 4603 (1985). 



